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Abstract 

It is shown that a satisfactory explanation of the ratio of the rates of the 





reactions 3 He(n~, nn) and 3 He(7r~ ,np) for stopped pions is obtained once the 
effect of the short range two-nucleon components of the axial charge operator 
for the nuclear system is taken into account. By employing realistic models 
for the nucleon-nucleon interaction in the construction of these components 
of the axial charge operator, the predicted ratios agree with the empirical 
value to within 10-20%. 
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1. Introduction 

A natural ansatz for the description of nuclear pion absorption and pro- 
duction reactions in the near threshold regime, where the reactions mainly 
involve S'-wave pions, is to express the reaction amplitudes in terms of the 
axial charge density operator of the nuclear system [1]. This operator is 
formed of a single nucleon term and a two-nucleon "exchange current" term. 
The latter separates into a long range pion exchange term, first derived in [2], 
and a short range term, the most important part of which can be determined 
directly from the nucleon- nucleon interaction [3]. Direct evidence for the 
short range component was recently found by the quantitatively successful 
explanation of the reaction pp ppir near threshold [4] , to which the main 
pion exchange term does not contribute [1] . This important role of the short 
range part of the axial exchange current operator has been substantiated in 
subsequent work [5,6]. 

In this paper we apply the ansatz of describing nuclear S'-wave pion ab- 
sorption by the axial charge operator to the case of absorption of stopped 
negative pions on 3 He. This reaction, for which the isospin structure of 
the nuclear states involved is richer than that in the case of the reaction 
pp — > ppn°, is less selective and therefore the rate also has an important long 
range pion exchange contribution. Yet if only this rescattering contribution 
is taken into account in addition to the single nucleon operator an overpre- 
diction by more than an order of magnitude of the empirical ratio [7] of the 
rates for production of nn and np pairs results. The point of this paper is 
to show that this problem is resolved by the effect of the short range con- 
tributions to the axial charge operator, which were recently found in Ref.[l] 
to provide the explanation for the cross section for the reaction pp — > ppn° 
near threshold. We show that when realistic models for the nucleon-nucleon 
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interaction are employed in the construction of the short range components 
of the axial exchange charge operator, a very satisfactory account of the re- 
cent data [7] on the predicted ratio of the two-nucleon absorption processes 
3 He(7i~ ,nn)p and 3 He(ii~ ,np)n is obtained. The short range corrections 
reduce the predicted value by the required large factor so that the final pre- 
dicted ratio exceeds the empirical value 6.3 ± 1.1 [7] by only 10-20%, the 
range of variation arising from the differences between the potential models 
considered. 

In section 2 we review the derivation of the S'-wave pion absorption op- 
erator, as it appears when constructed from the axial charge operator. The 
formalism for calculating the absorption rate for stopped pions by a pair of 
nucleons in a given eigenchannel is then presented in section 3. The numerical 
results and discussions are given in section 4. 

2. The S-wacve pion absorption operator 

We shall take the effective Lagrangian for the interaction between low 
energy pions and nuclei to have the form 

C(x) = -l-Apix) • &>4>(x), (1) 

J TV 

where 4> is the pion field, f n is the pion decay constant and the axial 
current density of the nuclear system. This Lagrangian density is the direct 
generalization of Weinberg's effective Lagrangian for the pion nucleon system 
[8] . In the case of a single nucleon the axial current operator has the form 

y 

A^O) =ig A u(p')^ 5 -u(p), (2) 

and thus in this case Eq.(l) is equivalent to the usual pseudovector nNN cou- 
pling model, since by the Goldberger-Treiman relation the coefficient gA/^fn 
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can be expressed as 



9a _ Unn /„>. 

where f n NN is the pseudovector coupling constant. 

The interaction Eq.(l) implies that S-wave pion production and absorp- 
tion at threshold is determined by the axial charge operator, which in the 
case of a single nucleon is 

A (0) = -g A (T-v^, (4) 

where v is the nucleon velocity operator. 

The most important two-nucleon contribution to the axial charge operator 
is found in Ref. [3] to be the 7rp-exchange operator, which has the form 

^-^ w 4(^) (t1xt2)+(i " 2) - (5) 

Here the symbol (1 <-> 2) represents interchange of all the coordinates of the 
nucleon pair. We denote the fractions of the momentum k of the absorbed 
pion that are imparted to the two nucleons as k\ and fc 2 (k — ki + fc 2 ). 
In the limit m p — > infinite the Tip exchange operator Eq.(5) reduces to the 
pion exchange operator of Ref. [2]. As noted in Ref. [9], it in this limit also 
corresponds to the isospin asymmetric term (r 1 x r 2 ) in the usual form for the 
S'-wave pion rescattering operator calculated using the usual pseudovector 
7riViV coupling and the conventional effective Hamiltonian for S'-wave irN 
scattering [10,11]: 



8ni f nNN er 2 ■ k 2 , Ai 2 . A 2 , , x 2 -i , /-, n \ (a\ 
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Here co>2 is the energy of the exchanged pion and uj n that of the initial pion. 
The coefficients Ai^ are the following combinations of the two S-wave pion- 
nucleon scattering lengths a 13 : 

Ai = ~m n (ai + 2a 3 ), A 2 = \m i: {ai - a 3 ). (7) 
o o 

The 7rp component Eq.(5) of the axial exchange charge operator is obtained 
by choosing the following values for the coupling constants [9]: 

Ai = 0, A 2 = — ^-^0.04, (8) 
16tt/^ 

which agree well with the values extracted from nN phase shift analyses 
[12,13]. 

The short range components of the axial exchange charge operator that 
were found to be important for the cross section of the reaction pp — > ppn 
near threshold correspond to two-nucleon operators that arise when the ax- 
ial field excites intermediate negative energy components of the nucleon 
(nucleon- antinucleon "pairs"), which are then de-excited by the short range 
components of the nucleon-nucleon interaction (Fig.l) [1]. Although these 
contributions are most readily derived by employment of explicit boson ex- 
change models, the corresponding operators can be constructed directly from 
any existing nucleon-nucleon interaction model, if that is expressed in terms 
of Lorentz invariant spin-amplitudes as e.g. the Fermi invariants. This con- 
struction has been presented in detail in Refs.[3,14]. 

In the decomposition of the nucleon-nucleon interaction in terms of Fermi 
invariants the short range components are contained in the amplitudes SVTA 
(the P-amplitude contains the long range pion exchange term, which has to 
be dropped once the np exchange contribution above is included explicitly in 
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the pion rescattering contribution). Numerically the most important of these 
short range axial exchange charge operators are those that are associated with 
the scalar (S) and vector (V) invariants: 



m 



N 



+^s(k 2 )r 1 x tV • k 2 } + (1 <-> 2), 



A °(V) = -^-{K(fc 2 )r 1 + v v {k 2 )T 2 ][<T l ■ P 2 + V x <t 2 ■ fc 2 ] 



+-^(fe 2 )r 1 x rV 1 • fc 2 } + (1 <-> 2), 



(9) 



Here we have defined the nucleon momentum operators P^ 2 for the two 
nucleons as (p 12 +p' 12 )/2 respectively. 

The potential components vs and vy can be given a direct interpretation 
as scalar and vector meson exchange potentials. In the case of a single meson 
exchange model these would have the expressions 

*«> = -iMh' (10) 

Here mf and are the masses of the isospin and 1 scalar and vector 
mesons respectively, and and gy are the corresponding scalar- and vector 
meson-nucleon coupling constants. Below we shall use the notation a and S 
to denote the isospin and 1 scalar mesons respectively. 
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Using the expressions above in defining the effective Lagrangian Eq.(f), 
the transition T-matrix for the absorption of a pion with charge a and mo- 
mentum A; by a nuclear system can now be written as 

S ft = 5 fi - i2n5 iA \P f - Pi - k) < | A a (k) | > (12) 

where the axial charge density is formed as the sum of one- and two- nucleon 
contributions 

A a (k) = A^(k) + A^(k). (13) 

We shall here restrict the expressions to the limit k — > 0. The choice of 
the kinematic variables for a two- nucleon system are indicated in Fig.l. For 
the one-body operator A^\k) we obtain from Eq.(4) the following matrix 
element in the momentum representation 

< P'l I I Pi > = C(2x) 3 (-l)«^ 

x (14) 
m N 

Here the coefficient C has been defined as 

C = i"*W , 1 UnN . (15) 

In terms of the kinematical variables in Fig.l, the contribution of pion 
S- wave rescattering term Eq.(6) to two-body axial exchange charge operator 
A^(k) is given by the following matrix element 

8-7T 

x j h- Tl _ M'lyV x t*]„} S ■ g + (1 « 2). (16) 



7 



Here v n represents the pion exchange interaction without coupling constants: 

vM = ~r^/ 2 M> ( 17 ) 

where F n (q) represents the vertex form factor, which here has been taken to 
be the same at the S- and the P-wave vertices. 

From the short range components of the axial charge operators defined 
in Eqs. (9) we obtian the following explicit contributions for the scalar and 
vector exchanges mechanisms (a, 00 for isoscalar and 5, p for isovector) : 

<p[p' 2 \A^(k)\p 1 p 2 > = C(-irr 1 a ^v a (q) ( T 1 .P 1 + (1^2)(18) 

m N 

<p\p' 2 \A^\k)\p 1 p 2 > = Ci-lYrl^-vsiq)* 1 ■ P, + (1 <-> 2){19) 



m 



N 



<p' 1 p' 2 \A^\k)\p 1 p 2 > = Ci-lTrl-^vM) 



m 



N 



x {a 1 ■ P 2 + % -<j x x cr 2 • q} + (1 <-> 2), (20) 



<p[p' 2 \A^(k)\p l p 2 > = C{-\r- 2 -v p {q) 

m N 



X {rlicr 1 ■ P 2 + V x a 2 ■ q} 
- ^xrV'^+fl-S). (21) 
Here we have used the notation 

*<«> = i0k F ' (qf - (22) 

where Si = —1 for i — a, 5, Si — +1 for i = u,p , and the g^s are the 
corresponding meson-nucleon-nucleon coupling constants . In the case of 
the Bonn potential model [15] the form factors and coupling constants are 
determined by a fit to nucleon-nucleon scattering data (We here use the 
parameters values in Table A. 3 of Ref. [15]). Note that the terms due to the 
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tensor coupling of p and u vector meson were found [3] to be unimportant 
and hence are neglected here. 

For a general nucleon-nucleon potential, such as the Paris potential [16] 
used in this work, the derivation of axial charge operator can be done [3,14] 
by the expansion of the potential in terms of Fermi invariants. The re- 
sulting expressions are the same as above except that the radial functions 
v<r(o),vs(q),v u (q), and v p (q) are replaced respectively by v£(q),Vs(q),v v (q), 
and Vy(q). In this work we consider the Paris potential [16] and use the v^(q) 
and Vy(q) calculated in Ref.[14]. 

The above expressions are used in the numerical calculations of n absorp- 
tion on 3 He at threshold reported below. 

3. Transition rate for the reaction 3 He(ir, NN) 
The kinematics of the reaction 3 He(n, NN) in the 3 He rest frame is illus- 
trated in Fig. 2. As we are interested only in the ratio between the absorption 
rates, which is not expected to depend sensitively on the wavefunctions, we 
shall neglect the nuclear distortion effect on the pion wavefunction. For sim- 
plicity the 3 He wavefunction is taken to be of a simple harmonic oscillator 
form with the oscillator parameter value b = 1.36 fm. It then has the sep- 
arable form ~ $(P)05 ;T (p), where P and p are the total and the relative 
momentum of a pair of nucleons in 3 He respectively. With these simplifica- 
tions and the kinematic variables shown in Fig. 2, the total transition rate of 
the reaction 3 He(7r, NN) is then determined by the following quantity 

Pm T1 m T2 ;a = r^E E / dPd P 5{u„{k) - 6 - -f- - -£-) | $(P) | 

x I E E < SsMsm S3 | jrrij >< TrMTm T3 \ tm t > 

S,T M s ,Mt 

x < p l ra Sl ra T1 ,p 2 m S2 ra T2 \ A \ (f)^ ,Mr ,k,a >| 2 , 
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(23) 

where r = s = 1/2, j = 1/2 and t — m t — 1/2 are the spin and isospin 
quantum numbers of 3 He, and obviously m S3 = rrij — Ms,m T3 = m t — Mt- 
The isospin quantum numbers and the momenta of the ejected two nucleons 
are (m Tl ,m T2 ) and p 1 = P + p/2,p 2 = P — p/2. The momentum and the 
charge of the pion are respectively denoted k and a. The outgoing two- 
nucleon state is properly antisymmetrized. 

In the stopped pion limit k — > Eq.(23) can be written in the following 
partial-wave decomposed form 

P mT1 m T2 ;a = < T ' M t I TTm T1 m T2 >< TrM T m rz I tm t > 

-7,(3 

x < T'M' T I TlM T a > g(~f,(3), (24) 

where 

9(1,0) = \j dP \ $ ( p ) H< 1,P\M h >W ( 25 ) 
with p = (m^ — e — 3P 2 /4M A r) 1 / 2 . Here we have introduced abbreviated 
partial- wave notations (3 : [LS]JT for the initial two-nucleon bound state, 
and 7 : [V S']J'T' for the ejected two nucleons. The matrix element of A has 
contributions from one-body term and two-body terms defined in section 2 

< 7 ,p | A | 4>p > = < 7 ,p | | <Pp > 

+ ^< T ,pM (,) \h> (26) 

i 

where i — ir, a, 5, p, cu. The contribution from the one-body term is 

< 1,P' I AM | $/3 >= -5(27r) 3 ^^0 [L5]JT (p), (27) 

where 4>[ls]jt(p) is a general bound state wavefunction. For an S-wave 
harmonic oscillator wavefunction of 3 He, we have L = 0, J = S and only 
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[S,T] = [0, 1], [1,0] are permitted by the Pauli principle. The radial wave- 
function is 4>[o,s]s,t(p) = 4>{p) — Ne~ b2p2 . The angular momentum coupling 
in Eq.(27) is isolated in the coefficient B. We find that 



B 



*—\—LL'S'(-iy j - s ' + V5j,j5 



Im 

~2 (2tt) 9 / 2 



MM' 



X 
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V 


L 
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J 
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J 
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S' 


1 

> 



(28) 



The matrix element of the pion S'-wave rescattering term is 



x |Ai < ST' || Oi || ST > -f 
where 



A 2 < ST || Os II ST > 



(29) 



with 



4&(0 = / P 2 « ( ^ ) (P , ,P)0[L5]JT(P)P, 
BP L {l)=p' J p 2 dpF^{p\p)4, {LS \j T {p), 

fW(j/,p) = (2tt) J^P^dxv^p'-p), 



(30) 



(31) 



where x = p' ■ p. Defining similar partial-wave projected integrals for the 
i = a,5,p and u; mesons exchange contributions, we then get the following 
expressions 

fiU 1 



< 7 ,p'Mi CT) \h> 



B- 



2m 2 N 

x < S'T' || Oi || ST >, 



{aW(lo + s^(l)} 



< 7 ,P' |0/3> = 5 



gA 1 

2/. 2m 2 



r {Aa.(L')+S«(L)} 



x < ST || 2 || ST >, 
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(32) 



(33) 



>= B 



9a 



2f n 2m 2 N 

x{-(AP L (L') + BP L (L)) < S'T' || Ox || ST > 
+V2(A^ L (L') - BP L (L)) < S'T' || 4 || ST >}, 



(34) 



x {-(^W (L') + BW(L)) < S'T' || 2 || ST > 
+ v^(AW (L ; ) - fl{&(L))(< S'T' I) 4 || ST > 
+ < S'T' \ \0 3 \ \ ST >)}. (35) 

The reduced matrix elements in the above equations are defined by the fol- 
lowing spin-isospin operators 



1 = tV-tV, 

2 = rV 1 - rV 2 , 



o ? 



V2 



r'xTV + ff 2 ) 



4 = (t 1 + t 2 )^ CT 1 x ( t 2 . 

v2 

It is then easy to show that 

<S'T'\\0 1 \\ST> = 6ST ((-1)( 5+T ) - (-1)( S '+ T ')) 

T" T 1 1 J S" 5 1 
1/2 1/2 1/2 J | 1/2 1/2 1/2 



(36) 



,(37) 



< S'T' \\ 2 \\ ST > = QST ((-1)( 5+T ') - (-1)( 5 ' +T )) 

V T \ \ [ S' S 



x 



1/2 1/2 1/2 1 1/2 1/2 1/2 



,(38) 
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< S'T II On II ST> 



x < 



eVTssf ((-i) 5 + (-if 

T T 1 
1/2 1/2 1 
1/2 1/2 1 



> < 



S' S 1 
1/2 1/2 1/2 



<S'T' \\0 4 \\ST> = 6y/lSST ((-1) T + (-if) 

5' 5 1 



x < 



1/2 1/2 1 
1/2 1/2 1 



V T 1 

1/2 1/2 1/2 



(39) 



• (40) 



4. Results and Discussion 

It is interesting to first note that all transition matrix elements , Eq.(27) 
and Eqs.(29)-(35), depend on the same geometric factor B defined in Eq.(28). 
Evaluating this factor for an initial S'-wave (L = 0) NN pair, it is easy to see 
that only the following transitions can occur in the reaction 3 He(ir, NN): 

fu = g^Po^So) 
f w = g( 3 P l , 3 S 1 ). 

Here the rate factor g is that defined in Eq.(25), and we use the standard 
spectroscopic notation 2S+1 Lj. The subindices on the f's indicate the total 
isospin T' and T of the final and initial two-nucleon systems. Evaluation of 
the Clebsch- Gordon coefficients in Eq.(24) leads to the following ratio of the 
absorption probabilities (Pm T1 m T2 ;a defined in Eq.(24)) for negative pions : 

3 He(n-,nn) = P-i/ 2 ,-i/ 2; -i 

5 He(7T-,np) P-l/ 2 ,l/2;-l 

= 1 + 6^. (41) 
in 
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The pion rescattering term A^\ Eq.(29), gives a very small contribution 
to /n because its main (A2) term depends on the isospin operator r 1 x r 2 
and therefore its matrix element between two T = 1 states vanishes. The 
remaining term gives only a small contribution because the coefficient A x ~ 0, 
as discussed in section 3. On the other hand, the pion rescattering term 
gives a large contribution to the transition between T = 1 and T = states 
( i.e. /10) and hence the prediction based on the sum of the single nucleon 
contribution and the pion rescattering term A^ is very large, as shown 
in the first row of Table 1 . 

In Table 1 we also show the effects caused by the short-range absorption 
mechanisms as calculated using the Bonn [15] and Paris [16] potential mod- 
els. By comparing rows 1-3, it is seen that the short range corrections due 
to the a and u exchange mechanisms (or more generally from the isospin 
independent scalar and vector exchange mechanisms) drastically reduce the 
predicted ratio. The contributions from isovector 5 and p exchanges are less 
important, but are significant in further reducing the predicted ratio towards 
the experimental value [7]. Although the contributions from the individual 
two-nucleon mechanisms differ between the two potential models, the net 
result obtained by taking all of them into account are quite similar and close 
to the experimental value. This situation is in agreement with that found 
for the reaction pp — > ppir° near threshold [1]. 

The results presented in Table 1 further strengthens the case for the 
importance of the short range pion absorption mechanisms [1]. It will be 
interesting to explore whether this model also is able to resolve several long- 
standing problems associated with nuclear pion absorption at higher energies. 
In the case of the two-nucleon system, all existing unitary 7riViV calculations 
[17] based on pion rescattering mechanisms ( including S and P irN partial 
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waves) fail to give satisfactory descriptions of the analyzing powers for the 
reactions n + d — > pp and NN — > itNN. Although the short-range absorption 
mechanisms that involve intermediate AApairs ( Fig.l) is in general weaker 
than the pion rescattering term, it can nevertheless have large effects on the 
polarization observables through interference with the large pion rescattering 
amplitude. This can be explored straightforwardly within the unitary nNN 
formalism developed by Lee and Matsuyama [18]. 

It will be also interesting to investigate the effect of the short-range ab- 
sorption mechanisms on 3 He(ir~ ,np) reaction in the A region. This reaction 
involves small contribution from the A excitation, since the dominant NA 
S'-wave intermediate state is excluded [17] for an initial 1 S'o pp pair in 3 He. 
The angular distributions of this reaction as well as the energy-dependence 
of the ratio of 3 He(ir + ,pp)/ 3 He(n~ ,np) have not been well understood [17]. 
An extension of the theoretical calculations, such that developed in Ref.[18], 
to include the short-range absorption mechanisms is expected to significantly 
change the theoretical predictions of these absorption observables for 3 He. 
A rigorous theoretical calculation of two-body absorption on 3 He is essential 
for extracting from the data the contribution from three-body absorption 
mechanisms which have been empirically established in Ref.[17]. 
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Table 1 



The ratios R = 3 He(rr~ ,nn)/ 3 He(ii~ ,np) as predicted using the Bonn 
[11] and Paris [12] potential models in the construction of the two-nucleon 
operators are compared with the experimental value [7]. The absorption 
mechanisms included in each calculation, defined in Eq.(27) and Eqs.(29)- 
(35) , are indicated in each row. 



Mechanism 


R-Bonn 


Rparis 


R-Exp. 


A W + A« 


533.42 


533.42 




,4(1) _|_ A^ +a ^ 


74.30 


208.12 






12.40 


16.81 




AW + A^ +a+U}+s ^ 


11.77 


11.88 




AW _|_ j\(+1T+a+U+6+p^ 


9.786 


8.0349 


6.3 ± 1.1 



Figure Captions 
Fig. 1: Axial exchange current operator arising from intermediate NN 
pairs that contributes to the NN <-> NNn reaction. 

Fig. 2: Kinematics of 3 He(ir, NN) reaction in the rest frame of 3 He. 
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